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Abstract 

Dynamics of planar domains with moving boundaries driven by the 
gradient of a scalar field that satisfies an elliptic PDE is studied. 
We consider the question: For which kind of PDEs the domains 
are algebraic, provided the field has singularities at a fixed point 
inside the domain? The construction reveals a direct connection with 
the theory of the Calogero-Moser systems related to finite reflection 
groups and their integrable deformations. 

1 Introduction, non-homogeneous porous medium 
flows 

Laplacian growth is a process that governs the dynamics of the boundary 
dQ = dtt(t) in the plane separating two disjoint, open regions (f2 and C\f2) 
in which harmonic (scalar) fields are defined. These may be interpreted as the 
pressure fields for incompressible fluids (porous-medium flows, Hele-Shaw flows, 
etc. See e.g. [11] and references therein). In the recent literature there have 
appeared several new formulations of the Laplacian growth related with the 
theory of integrable systems and random matrices, quantum Hall effect and 
Dirichlet boundary problem in two-dimensions (see e.g. [6], [7], [12] and refer- 
ences therein). 

In the present paper, we consider a new connection with the theory of quan- 
tum integrable systems. To be more precise, we study an integrable generaliza- 
tion of the Laplacian growth, when the boundary is driven by a field satisfying 
an elliptic PDE, that is not generally reduced to a Beltrami-Laplace equation 
("Non-Laplacian" growth). To be specific, we use the porous medium fluid 
dynamics interpretation. 

We find variable-coefficient elliptic PDEs for which the boundary dynamics 
can be described explicitly and the moving fluid occupies evolving algebraic do- 
mains (see below). These, turn out, to be PDEs of the Calogero-Moser type, 
related to finite reflection (Coxeter) groups as well as their integrable deforma- 
tions (see Conclusion) , that possibly complete the list of all second-order PDEs 
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connected with the algebraic domains (our main conjecture). 

In this section we set the problem in terms of the porous medium fluid 
dynamics. Formulation of the problem in terms of quadrature domains is given 
in the next section. 

Consider a flow of a liquid in a thin non-planar layer of non-homogeneous 
porous medium. The layer can be viewed as a 2-dimensional surface embedded in 
the three-dimensional Euclidean space. We let the layer curvature, permeability, 
porosity and thickness depend on the surface spatial coordinates x, y. We can 
choose x, y such that locally 

dl 2 = G(x,y)(dx 2 + dy 2 ) 

where dl is the surface length element. The surface area element is dT, = Gdxdy, 
and the volume of the liquid that can be absorbed in the range x + dx, y + dy 
equals 

dV — rjhdT, — rjhGdxdy, 

where rj = rj(x, y), h = h(x, y) are the medium porosity and the layer thickness, 
respectively. 

In the porous medium, the flow velocity v = (dx / dt , dy / dt) is proportional 
to the gradient V = (d/dx,d/dy) of the pressure P 

v = ^=VP, 

Vg 

where k = k(x, y) is the medium permeability. 

It is seen from the above that only the two combinations of variable coeffi- 
cients, namely 

vhG ' 73 

enter the flow equation of motion, and it is convenient to absorb h and G into 
definitions of the other coefficients. Therefore, without loss of generality, we 
can consider the flow in the plane parametrized by the complex coordinates 
z = x + iy,z = x — iy, choosing r\ and k to depend on z, z, while setting 
remaining coefficients to unity. The liquid volume conservation leads to the 
continuity equation 

(V • rru) = 0, (1) 
while the dynamical law of motion rewrites as 

v = -kVP. (2) 

We consider a situation where the liquid occupies a bounded, simply-connected 
open region of the plane, whose time evolution = Q(t) is induced by the 
flow. 
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At fixed time t, the pressure is constant along the boundary 

P(dSl(t)) = P (t). (3) 

Note, that the in the case of the simply-connected domains considered here, the 
dynamics is independent of Po{t), and for convenience the latter can be set to 
zero. 

The normal velocity of the boundary v n and that of the flow coincide at <9f2 
v n = n ■ v if z £ dQ. (4) 
The flow is singularity driven. For instance 

o k -«(*) ,_u „ \,\^( ^'(*) 



2k(zi, zi)?7(zi 



as z — > zi, 



(5) 

when a multipole source of order fc + 1 is located at z = zi. Equations (1) - 
(5) constitute a free boundary problem where the evolution of the boundary 
dil(t) is completely determined by the initial condition dfl(0) and strengths 
q{t) = q(t), fj,j(t),fij(t),j = 1, .., k as well as position z\ of the sources. 

2 Conservation Laws and Quadrature Domains 

From (1), (2), (5) it follows that the pressure satisfies the elliptic PDE 

VK,r ] VP = -Kq[S(x-x 1 )5(y-y 1 )]. (6) 
where q = q(t) is the differential operator of order k 



j=i \ / 



(7) 



Let (j>(z, z) be a time-independent function satisfying 

Vk7?V0 = 0, zeSJ (8) 

in whole Q, including z = z\. 

Let us now estimate the time derivatives of the following quantities 



M[<f>]= / rj(j)dxdy. 
Jn(t) 



Considering an infinitesimal variation of the fluid domain Q(t) — > Cl(t + dt), we 
get 

dM[<f>] 



dt 



v n r]4>dl, 
dn(t) 
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where dl is the boundary arc length. From (2), (3), (4) it follows 
dM[4>] 

ldfl(t) 



dt 



<j> (PKr)V<j> - <j)KT]VP) ■ ndl. 
Jdn(t) 



Applying the Stokes theorem and remembering that P and </> satisfy (6), (8), 
we get 

^ =*«(*)[*](*,*). (9) 

Note that mixed derivatives are absent q (c.f. (7)), for by (8), is expressed 
through first derivatives of <j>. q 3 - is the complex conjugate of qj, since both 
4>{z,z) and <j>(z,z) satisfy (8). 

It follows that M\<j>] is conserved for any solution of (8), a such that q[4>](zi,zi) 

0. 

The conservation laws have been first obtained for the homogeneous medium 
flows in [9] , the variable-coefficient generalization seems to be first presented in 
[5]. 

The flow in the homogeneous medium 

K = l, 7] = 1 (10) 

is the simplest example, where the conservation laws can be written down ex- 
plicitly [9]. In this example, any (anti)analytic in ft function satisfies (8) 

4>(z,z) = f(z) +g(z), for K = l, T} = 1, (11) 
where /, g are univalent in and the quantities 

/ (f{z){z- Zl ) k +g(z){z-z l ) k )dxdy 
Jn(t) 

are integrals of motion for the free-boundary flows driven by a multipole source 
of order k + 1 located at z = z\ in homogeneous medium. 
Returning to the general case, we integrate (9) getting 

M[0](t) = M[0](O)+7rQ^]( a! i,zi), 

where 

^^(tV^q + f^ + Q^j, . (12) 

Therefore M[(f>](t), and consequently a form of the domain, does not depend on 
the history of the sources and is a function of "multipole fluxes" 

Q= [ q(t')dt',, Q 3 = f q 3 (t')dt', Qj= f q 3 (t')dt', j = l..k 
Jo Jo Jo 
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injected by time t. 

Now consider the special case when M[0](O) = that describes the injection 
of the fluid to an initially empty medium. In such a case 

r)(z, z)(f>(z, z)dxdy = nQ[<j>(z k ,z k )], (13) 

Equation (13) is a generalization of quadrature identities (those expressing in- 
tegrals over 17 through evaluation of integrands and a finite number of their 
derivatives at a finite number of points inside Cl) appearing in the theory of har- 
monic functions [10] to the case of elliptic equations with variable coefficients. 
Special domains for which the quadrature identities hold are called quadrature 
domains in the theory of the harmonic functions. We extend this definition to 
solutions of an arbitrary elliptic PDE with regular in il coefficients. 

To construct quadrature domains (or equivalently the domains formed by 
fluid injected into initially empty porous medium), one needs an explicit form 
of a general solution of (8). Such explicit solutions are available for a class of the 
second order differential equations that are related to the Schrodinger operators 
of integrable systems on the plane. We, however, postpone construction of these 
solution to Section 4 and first present our main result in the next section. 



3 The Main Result 

Let us start with the simplest possible example, where the liquid is injected in 
the initially empty homogeneous porous medium through the single monopole 
source at z = z\ = x\ + iy\. By symmetry of the problem the solution is a 
circular disc of the radius r(t), centered at z = z\ 

\z- zi| < r(t). (14) 

The pressure satisfies 

AP = -nq(t)5{x - Xl )S(y - Vl ), 
where the source power and the total flux are 

«(*) = -^, Q = r ' 

respectively. 

The remarkable fact is that the variable-coefficient problem, with the con- 
stant porosity and the permeability that varies as an inverse square of one 
Cartesian coordinate ^ 

K=-o, T] = l (15) 
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admits the same circular solution (14) if the flow is driven by a combination of 
the same monopole source (of strength q — q(t)) and a dipole source of strength 
qi = —qQ/2x\, both located at the point z = z\ 

Indeed, it is not difficult to check that 

p = r ~k i^ 2xix +p 2 + r2 ) l °gp - r x(yx p 2 — - p 2 + x ( x - xi ) ~ ( 2xix + p 2 ) lo s( r ) 

(17) 

where p— \z — Zi\, satisfies (16), and the boundary conditions (2), (15), 

dt x 2 \dn J p=r 

as well as (6) holds at the disc boundary. 

Also the following quadrature identity holds 



L 



4 

7rr 



4>dxdy = 7rr 2 ^(zi,zi) , 
\z- Zl \<r 4xi \dx 

which is a simplest non-trivial generalization of the mean value theorem for 
harmonic functions to the case of regular in (14) solutions of the elliptic PDE 

The above example is a special case of the main result presented below. 

Let fl(t) be a domain, resulted from the injection of the fluxes Q,Qj,j — l..k 
into a homogeneous medium through a multi-pole source of order k + 1 located 
at z = z\. Then the same domain can be formed by an injection of a special 
combination of fluxes Q,Qj,j — l..k through a multipole source of order k — 
(k + l)(s(n + l) + l) located at z = z\ into an initially empty nonhomogeneous 
medium with permeability 

- Tv n>l>0, s>0, r?=l (18) 



(z s + z s ) 2n (z s - z s ) 

and constant porosity, where s,n,l are integers. 

In more details, the multipole fluxes of nonhomogeneous medium problem 
must be fixed functions of fluxes of its homogeneous medium counterpart 

Q = Q, Qj = Qj(Q, Q u Q h Q 1; ...,<3fc,zi,zi), j = l..(fc+l)(s(n+/)+l)-l. 
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For instance, in the above example of circular k = solution in a medium with 
permeability (15) (that is the special case n = 1, I = 0, s = 1 of (18)) 

Q 2 

fc = (0 + l)(l(l + 0) + l)-l = l, Q = Q, Qi = -f~- 



Note that (18) can be rewritten in the form 
1 

C(x,y) 



/• = -r^, C(x,y) = J] («■ z ) m °, (a ■ b) := Rc(afe), 



a£K+ 

where 1Z = {a} is a set of root vectors of a finite reflection (Coxeter) group on 
the plane (Dihedral group). 1Z is invariant under reflections z — > z — 2 a in 
a mirror (line), normal to any root vector a € 1Z, and 1Z+ denotes a positive 
subset of 1Z, containing a half of all root vectors. The multiplicities m a must 
be non-negative integers that are functions on the group orbits. If I = 0, the 
reflection group has one orbit and m a = n. Otherwise, the group has two 
orbits and multiplicities m a take values n and I on each orbit respectively. The 
permeability (18) is therefore an invariant of the group of symmetries of a regular 
4s-polygon (2s-polygon if I = 0) and its singular locus coincides with the union 
of mirrors. 

The pressure satisfies the elliptic PDE 

VC(x, y)- 2 VP = C(x, y)- X HC{x, y)- x P = -*q[8{x - x k )8(y - y k )] 



order(<7) = (k + 1) (deg(C(x, y)) + 1) - 1 = (k + 1) 1 + ^ m Q - 1, 

where H is the Schrodinger operator of the Calogero-Moser system related to 
the dihedral group [8], [4], [2] 

jj * \ - sm a (m a + 1) 

H = A- ) (a -a) — -. r „ — . 

^ (a- z) 2 

Also, regular in £1 solutions <fi(z,z) of the elliptic PDEs V((x,y)~ 2 V<f> = 
satisfy the quadrature identities 

/ <j)dxdy = QMzuZ!), order(g) = (k + 1) (deg(((x,y)) + 1) - 1 
J si 

that are generalization of the mean value theorem for harmonic functions to 
the case of solutions of the variable-coefficient elliptic PDEs in (generally) non- 
circular domains. 

Before proving the main result we need to construct a complete set of solu- 
tions to (8) for the medium with K,r] given by (18). 
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4 Nonhomogeneous porous medium flows and 
integrable systems related to the finite reflec- 
tion groups 

In this section we show how to obtain a general solution <f>(z, z) to (8), (18). It 
is instructive to consider the special case s = 1, 1 = of (18) 

«=^T> ( 19 ) 

when permeability depends on one Cartesian coordinate only. Construction of 
solutions in the general case (18) is conceptually similar. 

Consider the simplest non-trivial example n = 1 in (19) and start with 
factorizing the differential operator d 2 as 

8 2 _ ( d_ 

dx 2 \xdx)\ dx 



By the associativity of the differential operators 

x ik- 1 ) \xdxJVdx~~ 1 ) = \ x dx~~ l )\xdi} ^sl -1 ) 



Therefore, 



T 1 A = L 1 T 1) T 1= x^--1, L 1= x 2 V^V (20) 
ox x z 

and the elliptic equation (8) is nontrivially related to the Laplace equation when 
kt] = 1/x 2 . In general, the identity 

TA = LT 

that relates two differential operators (e.g. A and L above) through a differential 
operator T is called an intertwining identity and T is an intertwining operator. 

A differential operator that is related to A through the intertwining identity 
equals, modulo a gauge transformation, a Schrodingcr operators of an integrable 
system. Indeed, 

T(A- A) = (L- A)T 

and when v is an eigenfunction of A with the eigenvalue A, T[v] is an eigenfunc- 
tion of L with the same eigenvalue or zero. 

The factorization approach leading to the simplest nontrivial intertwining 
identity (20) can be now applied to L\ etc. By induction we get the intertwining 
identity for an arbitrary nonnegative integer n in (19) 

T„A = L„T„, L n =x 2n V-^-V, (21) 
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where 

_ n (d_ _ n \(d_ _ n-l\ ( d_ _ 1\ _ y, { &_ 

\dx x) \dx x J \dx x) dx 1 ' 

Any solution to (8), (19) in Q can be represented in the form 

= T n [/], A/ = 0, zefl. (23) 

Let us show this for n = 1, where T\ is given by (20). Introduce / satisfying 
the following equation 

ox 

where Li[<f>] = in Q and <p = for z £ ft. It is not difficult to see that 



4>{x' + iy, x' - iy) ^ , 
(x' 



/(*,*)=*/ yV VAa yJ dx' + xF(y), 



where F(y) is an arbitrary regular function of y. Therefore, for any <fr regular in 
Vl (which is the case) there exist a regular in Q function /, such that Ti[f] = cj). 
From the intertwining identity (20) it follows that T\Af = and Af € Ker(Ti), 
if z e Q i.e. 

A[f]=xC(y), z = x + iyen 
where C(y) is an arbitrary function of y. It follows that 

A[f + xF(y)]=0, f!M = C(y). 

Since / is defined modulo xF(y), we can set F — and any solution of L\ [</>] = 
can be represented as <f> = Ti[f], where A[/] = 0. Similar proof applies to the 
arbitrary n case (23). 

In the general case (18) 

T n ,i ;s A = L n i. s T n i. s , L n i ;s = (z s + z s ) 2n (z s - z s ) 21 V — — — — — — —57V, 

(z s + z s ) zn (z s — z s ) u 

where the intertwining operator can be expressed in the form of a Wronskian 
[2] 



T M „'("+') W \- s[n ( i)> sin (^), -, sin(0 n ), /] 

(24) 

with 

e f fc( S 0+f), fc = l,2,..n-i 

p ' fc \ (2fc + /-n)(sfl+f) ) n-l<k<n ■ 
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It is important that the intertwining operator T n j- S is a homogeneous differential 
polynomial in z,z,d z ,dg. This fact allows one to construct the quadrature 
domains for solutions <p of L n ^ s [(f)] = 0. 

Note that in general, there exist several independent operators intertwining 
A and another second order differential operator (forming a linear space of inter- 
twining operators). For instance, T n in (21) that intertwines A with L n = L„ ; o ; i 
is not a special case T n q ; i of (24). However, any nonzero linear combination 
of them can be used to obtain solutions to the corresponding elliptic equations 

£n,O;l[0] = 0. 

5 Proof of the main result 

The zero- initial condition solution of the homogeneous-porous medium flow that 
is driven by a multipole source of order fc + 1 is described by the polynomial 
conformal map of degree fc + 1 from the unit disc in the parametric |w|-plane 
into the fluid region SI 

fc 

z{w) = zi + rw + ^2uiW i+1 , \w\<l (25) 
i=i 

The map is analytic in \w\ < 1, and the unit circle \w\ = 1 is mapped to the 
boundary dfl. As shown above such regions are also quadrature domains. They 
are special (polynomial) cases of (rational) algebraic domains [11]. The map 
coefficients r,u i7 i — 1, ..fc are functions of Q,Qi,Q i7 i — l..k. The "conformal 
radius" r can be chosen to be real. 

The main idea of the proof is to show that the quadrature identity (13) for 
solutions of (8), (18) holds in domains defined by (25). 

Consider illustrative examples of the problem (19) with permeability chang- 
ing in one direction. As shown in the previous Section any solution of the elliptic 
equation L n [<f>] = can be represented as 

T n [f(z)+g(z)}, (26) 

where f(z),g(z) are holomorphic and anti-holomorphic respectively. 
According to (13) we have to show that 

/ T n [f(z)]dxdy = nQT n [f(z)] z = Zl 
Jn 

holds for any analytic in SI function f(z) when SI is defined by the conformal 
map (25). Using the Green theorem and taking (22) into account we rewrite 
the last equation as 
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Since w = l/w if \w\ = 1, z(w) — z(l/w) along the boundary, and we can 
rewrite the left-hand side of the last equation as 

Y .L I ( ( Z (w)+z(l/w)y +1 ( [f( z ( w ))} \ ^Idw 

Since z(w) is analytic in |io| < 1, z(l/w) has poles only at w = 0, and the above 
integral is a pure sum of residues 

(fe+2)(n+l)-2 . • x 

E ^ . (») 



J=0 



where y,,j = 0..(fc+2)(n+l)— 2 are functions of the parameters zi,zi 7 r,Uj 7 Uj,j = 
l..k of the conformal map (25). Equating it with the right-hand side of (27), 
we see that Qk must be (differential operators) of order (k + l)(n + 1) — 1 and 



(fc+2)(n+l)-2 . ■ . 

QT„[/(z)] z=Zl = ^ C7 3 ; , (29) 



where C/j, j = 0..(fc + 2)(n + 1) — 2 are linear functions of Q, Qj,Qj,j = l..(k + 
l)(n+ 1) — 1. Therefore, the quadrature identity (27) is satisfied if the following 
system of2(fc + 2)(n + l) — 1 linear equations 

Vj~Uj=0, Vj-Uj = 0, j = 0..(fc + 2)(n + l)-2 (30) 

for 2(fc + l)(n + 1) unknowns 

Q,Q,Qj,Qj, j = l..(fe + l)(n+l)-l (31) 

has solutions. 

Note that the condition Q = Q is satisfied automatically, since, as easily 
seen from (27), (25), (12), the j = subset of (30) 

V - U = 0, V o -U o = 

are equations for Q,Q that have the same form for any n > in (22). They 
have real solution 

k 

Q = Q = r 2 +Y(i + l ) u i u i 

i=l 

The number of equations in (30) exceeds the number of unknowns (31) by 2n and 
the system of equations (30) is overdermined for the non-homogeneous medium 
problem n > 0. 
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For instance, for the circular domain k = 

z(w) = Z\ + rw 

in a medium with permeability l/x 2 (that has been considered in Section 3), 
system (30) consists of six equations 

Q - r 2 = 0, Q - r 2 = 0, 

2Qa;i - Qi + Qi - 2xir 2 = 0, 2Qa;i - Qi + Qi - 2x x r 2 = 0, 
4QiXi - r 4 = 0, 4QiXi - r 4 = 

for four unknowns Q,Q,Qi,Qi- It has the following solution (already demon- 
strated in Section 3) 

Q = Q = r 2 , Qi = Qi=r i /Ax 1 

Returning to the general case, we are going to show that not all equations in 
(30) are independent, the system is compatible and has a unique solution, which 
proves our main result presented in Section 3. 

To prove the compatibility, we introduce the basis 

<j> j (z,z) = c j T n [(z-z 1 y+ n ], j (z,z) = c j T n [(z-z 1 y+ n ], j = 0,1,2,... (32) 

where Cj = + n)!x™, in the space of solutions 4>(z, z) of L n [<p] — that are 
regular in neighborhood of z = z\, and then show that the quadrature identity 
(13) holds for any element of this basis. 

Indeed, according to (22), (32) is a set of solutions of L n [(j)] = that contin- 
uously tends to the basis of functions analytic in a neighborhood of z = z\ 

4>j{z,z) -> (z - zx) , j = 0,1,2,..., x x -> oo 

as the position of the source z\ = xi + iyi goes to infinity. On the other 
hand, L n [(f>(z, z)] = 0, z e O continuously tends to the Laplace equation when 
region is moved to infinity. Therefore, set (32) is homotopically equivalent to 
(z — z\) J , (z — zi) J , j = 0, 1, 2, .. under a continuous deformation of the Laplace 
operator and thus contains a basis of solutions of L n [(/>] = that are regular in 
a neighborhood of z = z\ . 

It then follows from (28) and (29) that the quadrature identity holds if 

^i-^icW) = °' ^^^-W) = °' ^=Mfc+2)(n+l)-2 

V / z—z\ V / z—z\ 

(33) 

for 

f(z)e{(z-z l y iJ >n} 

Since j > n, the left hand sides of the first n equations in (33) vanish, and 
there remains 2(k + l)(n + 1) independent equations and the equal number of 
unknowns (31). 
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We now have to prove the compatibility of remaining equations. (30) is 
a non-homogenous system of linear equations for unknowns Q,Q,Qj,Qj,j — 
l..(fc + l)(n + l) — 1, that fixes dependence of these unknowns on parameters 
z\, z\, r, Ui, u^, Mi, Wj.. The system is compatible if its homogenous part 
does not have nontrivial solutions. Let us suppose that it does. Remind that 
the homogenous part of the system has been obtained by action of the operator 
Q to an arbitrary solution of L n [<p] — at point z = Z\. So, if the homogeneous 
part of the system had nontrivial solutions, then operator Q would annihilate 
any solution 4> of L n [(j)] = at z = z\, i.e. 

QT n [f(z)] = at z = z u (34) 

where f(z) is any analytic in f2 function. If the above were true, then changing 
continuously the position z = Z\ of the source, we could construct such operator 
Q, with coefficients depending on z, that QT n [f(z)] = in some region of the 
plane for an arbitrary f(z). But this is evidently impossible, since the highest 
symbols of Q and T n contain pure derivatives in z, so does their composition. 

Therefore (30) has a unique solution and quadrature identity holds in the 
polynomial algebraic domains, that leads to our main result (see Section 3). 

Similar result can be analogously proved for the general system with perme- 
ability given by (18). 

6 Generalization to Higher Dimensions 

Our study can be extended to PDEs in more than two dimensions, since the 
derivation of the conservation laws in Section 2 is not restricted to the two- 
dimensional flows. 

For instance, if we take solution 0(£) of V^ 1 " 2 V0=O, where £ :— (£i, ..,£(;) 
and V := ■•, , then the following quadrature identity holds 

i|€-€'l<r V (a + 2)4i OKi / 5= ^ 

in the (i-dimensional ball of radius r with center at £ = where v d = L^_^n <r 
is volume of the ball. Since the Hadamard expansion of the fundamental so- 
lution of the Calogero-Moser type operators truncates (see e.g. [3], [4]), the 
pressure distribution in the corresponding free-boundary flow can be written 
down explicitly for any d (e.g. we have used technique of [3] for derivation of 
(17))- 

7 Conclusion, Main Conjecture 

In this article we have found examples of non-homogeneous porous medium 
flows, driven by a multipole source at a fixed point, whose boundaries obey 
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the same dynamics as those of the homogeneous-medium flows also driven by 
a multipole source located at the same point. Namely, the medium permeabil- 
ity is a homogeneous rational function of x, y, and an invariant of the dihedral 
group. The multipole fluxes of the non-homogeneous medium problem must be 
fixed functions of those of the homogeneous medium one. The related variable- 
coefficient elliptic PDEs for the pressure distribution are of the integer multi- 
plicity Calogero-Moser type and the quadrature identities for solutions of such 
equations hold in polynomial algebraic domains. 

The above result has been derived using technique of intertwining operators. 
The fact that coefficients of these operators are polynomials in z,z is essential 
for algebraicity of the corresponding quadrature domains. This result can be 
extended in the following way 

According to works on algebraic integrability [4] , it is likely that the second- 
order elliptic operators related to the Coxeter root systems as well as their 
special deformations exhaust all possible operators that can be related to the 
Laplace operators through polynomial intertwining operators. These are ellip- 
tic operators for which a polynomial Baker-Akhieser function exists. We call 
all such operators the Algebraic Calogero-Moser operators. For instance, the 
algebraic Calogero-Moser equations [4] related to the porous medium problems 
(8) with the permeability 

k = x~ 2m ((2m + l)y 2 - x 2 )~ 2 , r] = 1, 

where m = 2, 3, 4,..., are simplest examples that extend the main result of the 
present paper to no n- Coxeter arrangements of mirrors. 

More generally, the problem of classification of all such algebraic Calogero- 
Moser systems in two dimensions is as follows [1], [2], [4]: Find a strictly increas- 
ing sequence of integer positive numbers < ki < fc 2 < . . < k n and a sequence 
Ldi, ..oj n of complex parameters ("phases"), such that the ratio of Wronskians 



nkn lF[ sin(fl 1 ),..sin(fl„_ 1 ),sin(fl w )] 
W[sin(0i),..sin(0 n _i)] 



C(x, y) = - L - T ^r'n \ n . e^W + uj, z = pe i ' 



is a polynomial in x, y. The corresponding intertwining operator is of the nth 
order and a polynomial in x, y that can be rewritten in the form of a Wronskian 



T[f] = P 



kn W[zm(0{),...zm{Q n ),f] 



W[sin(0i),...sin(0 n _i)] 



with (24) being a special case. The results of the present paper (i.e. algebraicity 
of domains for PDEs with n = l/£(x,y) 2 ) also hold for the deformed systems. 

In conclusion, one may pose the classification problem: Find complete list of 
all PDEs whose solutions satisfy quadrature identities in polynomial algebraic 
domains. In view of the above it is reasonable to expect that it has the following 
solution: 
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Conjecture: The algebraic Calogero-Moser equations exhaust, up to a gauge 
equivalence, all possible second order elliptic PDEs whose solutions satisfy quadra- 
ture identities in polynomial algebraic domains 

As mentioned in Section 6, our study can be extended to PDEs in more 
than two dimensions. Although, we cannot use conformal maps to parametrize 
algebraic domains in higher dimensions, we can still pose the classification prob- 
lem for domains with spherical boundaries, looking for a complete list of elliptic 
equations for which generalized mean value theorem holds. More precisely, one 
can look for equations in c?-dimcnsions whose solutions possess the following 
property: Integral of an arbitrary solution taken over an arbitrary d-dimensional 
ball, equals a linear combination (with coefficients depending on the ball radius 
and position) of the value of the solution and those of a finite number of its 
derivatives at the ball center. 

Extending our conjecture to higher dimensions we may expect that algebraic 
Calogero-Moser equations complete the above list. Note, that our classification 
problem in d > 2 seems to be equivalent to a restricted classical Hadamard prob- 
lem for irreducible Huygens' operators that non-trivially depend on more than 
two variables. As in our case, all known examples of such Huygens' operators 
are related to algebraic Calogero-Moser systems [4] . 
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